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Abstract 

Relativistic heavy ion collisions can generate metastable domains in which the chiral 
condensate is disoriented. Nucleus-sized domains can yield measurable fluctuations in 
the number of neutral and charged pions. We propose a scenario in which domains 
are 'annealed' by a dynamically evolving effective potential in the heavy ion system. 
Domains of sizes exceeding 3 fm are possible in this scenario. 



The order of the chiral restoration phase transition in QCD for reahstic values of the 
up, down and strange quark masses is currently unknown [|l], 0]. If the transition is nearly 
second order, Rajagopal and Wilczek have speculated that the nonequilibrium dynamics 
of a heavy ion collision can generate metastable domains in which the chiral condensate is 
disoriented. Such domains can provide a coherent source of pions that can exhibit novel 
"centauro-like" fluctuations of neutral and charged pions [|], |^, 0- However, the ability of 
the coming ion-ion experiments at the RHIC and LHC colliders to resolve such fluctuations 
critically depends on the size and energy content of these domains. 

Numerical simulations [§] show that the zero-temperature linear sigma model advocated 
in leads to domains that are roughly pion sized, too small to be resolved in experiments. 
Nevertheless, we illustrate in this Letter that finite temperature effects in a heavy ion collision 
can enhance the size of domains relative to the zero-temperature estimates. 

Rajagopal and Wilczek propose a "quench" scenario in which the condensate is initially 
chirally symmetric as appropriate at high temperature, but its evolution is taken to follow 
classical equations of motion in the absence of a heat bath.0 For this quench scenario to be 
realistic, the expansion rate of the heat bath needed to create the symmetric initial state 
must greatly exceed the rate at which the mean field evolves. The opposite is more likely the 
case. In and p], the a and vr fields in the 0(4) linear sigma model are taken to characterize 
the dynamics of the quark antiquark condensate in two- flavor QCD. These fields evolve over 
the time scale ~ ~ 1/3 fm ||^, where m„ = 600 MeV is the zero-temperature a mass. 
On the other hand, model calculations |^ indicate that chiral symmetry is broken very late 
in RHIC-energy collisions, perhaps at proper times Tc ~ 5 — 20 fm, so that the expansion 
time scale ~ Tc is much larger in comparison. 

We present an alternative "annealing" scenario in which the sigma-model condensate 
evolves semiclassically in the presence of a bath of quasiparticle excitations. Fluctuations 
induced by the quasiparticles create an effective potential in which a chirally symmetric 
state is initially stable (modulo finite quark-mass effects). As the system expands and 
rarefies, the effective potential develops a "wine-bottle" shape and slowly changes towards 
its free-space depth. The condensate then evolves in this changing potential. 

Transient domains of disoriented $ = (a, vr) condensate can develop because the $ ^ 
initial state is unstable in the wine bottle potential. Consider first a quench scenario in 
which the potential, 

Vo{^) = A(<l>2 - v^)/A - Ha, (1) 

is fixed. The system "rolls down" from the unstable local maximum of Vo($) towards the 
nearly stable values with |$| = v, since the symmetry breaking term —Ha is relatively 
small. This process is analogous to spinodal decomposition in condensed matter physics 
n], 0. Field configurations with vr 7^ develop during the roll-down period. The field 



will eventually settle into stable oscillations about the unique vacuum (/^, 0) for H ^ 0, but 

^ This procedure is not the standard quench in condensed matter physics, in which the system remains in 
contact with a heat bath of zero temperature. We continue to use the term "quench" here, however, because 
a T = heat bath would also effectively imply decoupling from the system since the interactions of pions 
effectively disappear at zero momentum due to the approximate chiral symmetry. 
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oscillations will continue until interactions eventually damp the motion. In the heavy-ion 
system a domain can radiate pions preferentially according to its isospin content. 



Domain growth in the annealing scenario with an evolving effective potential differs from 
the quench as follows. Initially, is nearly flat for a ~ vf as shown in fig. 1. Therefore, 
the roll-down time scale can be very large at first, or infinite if the corresponding equilibrium 
transition is truly second order (as is certainly the case for H = 0). Only as the potential 
approaches its free space shape does the roll-down become rapid. In both quenching and 
annealing, domains can grow as long as a is substantially different from /^r. In the annealing 
scenario, however, that time scale is limited by the time needed for Ves to approach Vq. That 
time in turn depends on the global dynamics of the nuclear collision. 

We describe the evolution of the condensate in the context of a self consistent Hartree- 
like semiclassical approach. There, the condensate evolves as a classical field in the presence 
of a fluctuating quasiparticle bath, see e.g. |]TB[. The condensate fields $ = (cr, tt) obey 



relativistic Ginzberg-Landau |]T4[ equations: 

{dyOt^ - V')$ = X{v^- 3{6^) - {6<^l) - $ - Hn^, (2) 

to quadratic order in the fluctuations. The symmetry breaking field H is fixed in the sigma 
direction, n„ = (1,0). Both the condensate $ and the fluctuation fields 5$ can instan- 
taneously have any orientation in the 0(4) internal space. When the T = condensate 
^ = (/tdO) is realised, the fluctuations 6^\\ along $ can be identified with the sigma field 
with the familiar tree-level mass rria- = (2A/^ + m^)^/^. Similarly, the three transverse modes 
are pion-like with the T = mass = (H/ f^^Y^'^. 

The Hartree correction 3(5$y) + accounts for the interactions of sigma-like and 

pion-like quasiparticles with the condensate. To estimate this correction, we must understand 
the nonequilibrium dynamics of the quasiparticles in the heavy ion system — no trivial 
task! To simplify our exploratory treatment, we assume that the quasiparticles are localized 
excitations of momentum p and energy EZ'^ = [p^ + ('^eff^)^]^''^) where m"^ are the effective 
masses of the sigma- and pion-like quasiparticles. We can then write the Hartree term in 
terms of semiclassical phase space distributions fa^^T,{p,r,t): 



3(5$J) + = 3 / dT,^ + {N-1)J dvj-^ (3) 



where dVp = d^p/{2TiY and = 4. Note that for quasiparticles near local thermal equi- 
librium, the /a(p) are Bose distributions so that 3(5<l>|) + (5$^) ~ T^/2 for T ^ m"^ and 
m-gg. The equations of motion (|]) with (0) then correspond to the more familiar effective 
potential shown in fig. 1. The critical temperature at which the wine bottle shape of Ks 



disappears is then Tc = V^fn ~ 132 MeV. [Lattice simulations suggest that the condensate 
in equilibrium QCD substantially changes at temperatures that are somewhat higher than 
this 0(4) model estimate 



In general, the quasiparticle distributions will satisfy Boltzmann-Vlasov kinetic equations 



in which m'^g are self consistent functions of the field (see, e.g. [0). Here we will neglect 
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the following: (z) the effect of collisions among quasiparticles, and (ii) the dynamical effect of 
the quasiparticle masses. It is reasonable to neglect collisions since, by the time T < T^, the 
system is sufficiently dilute that the collision frequency is small compared to the expansion 
rate [0. On the other hand, assumption {ii) is certainly wrong, since m'^Q must eventually 
change from ~ 171.,^ at Tc to its free-space value ^ 600 MeV. f\ We emphasize that these 
approximations, while crude, allow for the most rapid "quench" consistent with causality, 
since the quasiparticles then stream away, unimpeded by the collisions or the gain in mass. 

To illustrate how the fields and quasiparticles evolve with time, we restrict our attention 
for the moment to a single domain in the interior of the collision volume. An expanding 
system created in a central Au+Au collision at RHIC will reach temperatures below ~ 
132 MeV only after a time ~ 5 — 20 fm comparable to the Au transverse radius Ra ~ 7 fm. 
We take the quasiparticle flow to be roughly 3-dimensional and homologous, allowing for the 
fastest (causal) flow. We assume that at this late time in the evolution of the reaction, the 
condensate is invariant under radial boosts, so that it is a function only of the proper time 
r = (t^ — r^)^/^ for radial distances r < Ra and late times t > Ra- The left side of (0) is 
then 

(4) 

Note that for the one-dimensional Bjorken expansion valid for r, t < Ra, one would replace 
r by (t^ — z'^Y^'^ and the "3" multiplying the flrst derivative by "1". However, we stress that 
chiral symmetry is unbroken within that space-time region in nuclear collisions. 

Similarly, the quasiparticle distribution functions r, t) for a = a, n are functions 
only of p' = {pr — Eft)/T and r. Approximation (z) implies that f{p) satisfles a coUisionless 
Boltzmann equation (see and refs. therein) describing the free-streaming evolution of the 
phase-space distribution. We obtain fa{p,r,t) = /^{p't/tc), where the initial distribution 
faip) the time Tc is assumed to have the thermal equilibrium form [exp{Ep /Tc) — at the 
temperature = \/2/,r- Neglecting quasiparticle collisions, as done here, this form applies 
even if effective masses are functions of the condensate $(t). This solution also applies in 
quasiparticle-collision dominated regime where the expansion is adiabatic, provided that the 
mgfj can be neglected. The Hartree term (^) is then a function of r alone. With (ii), we 

(3(5$5) + (5$i))(r) ^ {3{6^) + (5$i))e(r,/r)^ (5) 

an approximation that is adequate provided that we focus only on the evolution of the 
condensate. 

We determine the average fleld in a single domain $ = (cr, tt) as a function of r by 
combining (0), (H) and (Bf) to flnd 



3d 



+ -_ \ ^ = Xv 



r ^2 



e I — 



$_A|$|2$-i7n,. (6) 



^By neglecting the time evolution of masses, we avoid the complicated treatment in that ensures that 
all quasiparticle modes have effective masses satisfying > 0. The present approach is adequate, since we 
focus on the behavior of the condensate, and not the quasiparticles. 
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The standard parameters A = 20, v = 87.4 MeV, and H = (119 MeV)^ are consistent with 
the values of the pion decay constant f^^ = 92.5 MeV and the zero temperature meson masses 
= 140 MeV and = 600 MeV. The anneahng parameter, 

e^.;-2(3(5<|.J),+ (5<|.i),), (7) 



determines the strength of the Hartree term. A value e = 1 applies for m'^g{Tc) -C T^. For 
m'^Q ~ 140 MeV at T^. ^ 132 MeV, we find e ~ 0.5. The value e = 1 applies if chiral 
restoration in an equilibrium system occurs as a true second order phase transition. On the 
other hand, a consistent treatment at the Hartree (tadpole) level for the 0(4) model with 
H ^ suggests m'^Q > and a gradual, continuous chiral restoration, rather than a phase 
transition. As mentioned earlier, the case of QCD applying to the real world is currently 
ambiguous. 

In Fig. 2a we show numerical solutions of (^) for Tc = 10 fm, tt^{tc) = 48 MeV, dir-^/ dT\r^ = 
95 MeV/fm, while in Fig. 2b we show the case vr'^(Tc) = 5 MeV, dn-^ / dT\r^ = 0.0 The dark 
lines and light lines correspond to calculations for e = 1 and 0.5, respectively. In fig 2a. we 
see that the time r/ needed for cr to begin to oscillate about a = f-„ is roughly 2.5 Tc for 
both values of e. As expected, tj is essentially the time needed for V^fj to fall within roughly 
20% of the T = potential. For e = 0.5, the solution for the more quiescent initial condition 
shows a similar structure, but reaches a ~ sooner, at rj ^ 1.3 r^. 

To study the size of domains, we now shift our focus from the average $(t) in a single 
domain to the spatial variation of the condensate due to fluctuations in the initial fields. 
Domains arise because a spinodal instability enables small fluctuations in the initial field 
configuration to grow [|12|. Neglecting the expansion of the system for the moment, we see 
that small fiuctuations about $ = in (^) "run away" as e™ where fl^ = — k"^, for a 
wave number k < r^"^. The k = mode grows the fastest, with a time scale 

T.,= [\{v'-3{S^l)-{S^l)]''^\ (8) 

which is a function of time r itself. In the heavy-ion system domains grow due to the 
macroscopic expansion as well as the microscopic instability of the infrared modes of the 
coherent chiral condensate. We use (0) to find that small fiuctuations follow 

= (9) 

for initial conditions specified at r = r^. The factor (tc/t)^ in front of fc^ in (y) describes the 
homologous expansion of domains; an analogous term enters in descriptions of the smoothing 
out of fiuctuations on a local scale due to the expansion of scales in infiationary cosmology 





In principle the initial conditions are determined by microscopic fluctuations. The initial spread in d^/dr 
and $ can be quite large compared to v. We estimate {(d^/dr)'^) = J2al d^p^pfc(p) ~ (95 MeV/fm)^ 
while ($2) ^J2al dTp{E^r^f^{p) « (48MeV)2. 
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To obtain a rough estimate of the growth of the domain size in the comoving frame, we 
follow Boyanovsky et al. 



ri|] and calculate the contribution of the growing fluctuations to 



the correlation function 



(7r°(f,t)7r0(0,t)) 



2^ 



ik-f 



(10) 



where (...) represents an average over the fluctuating initial conditions. For definiteness, we 
choose the vr" direction in isospin, but all cartesian pion directions are equivalent by isospin 
invariance. The initial fluctuation spectrum {ir^ik, rc)7r°(— fc, r^)) is thermal to a high degree 
of accuracy, yielding spatial correlations over characteristic distances of a thermal wavelength 
{7rT)-\ i.e. (7r°(f, t)7r°(0, t)) ~ exp{-TrT\f\). For {vrT}-^ < < Inpir^rfdr, we can 
then evaluate (H) using a WKB approximation and integrate (p!0D applying the saddle-point 
method. We find 



(7r°(r,t)7r°(0,t))^(7r°(0,t)')exp -r^ 



(rc/r)Vsp(r)d7 



-1 



This result suggests that domain size i?^ increases roughly as 



Roit) ~ 8 / (r,/r)Vsp(r)rfr 
Jo 





'sp 'c 



COS ^ f A/ey ) - COS ^(Ve) 



(12) 



where rP^ 



(Ai;2)-i/2 ~ 0.5 fm. 



Domains can grow as long as the system is unstable, i.e. until the a field begins to 
oscillate about the value a ^ f^^. For a nearly critical system with e = 1, we estimate this 
time to be Tf — Tc ^ 1.5 Tc, see figs. 2. Such values imply that the domain size can reach 

^ 7 fm — very large indeed! We obtain a more conservative estimate i?^ ~ 3 — 4 fm by 
taking e = 0.5 and supposing that domain growth stops after a reaches f^^ for the first time 
at Tf — Tc ^ (0.3 — 0.5)rc. Note that the numerical effect of the inflation factor {tc/t)"^ on 



(|T2|) is insignificant for the parameter values considered. Admittedly, numerical simulations 
for the annealing scenario are needed to provide more concrete estimates. Such simulations 
However, we point out that similar estimates worked surprisingly well 
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are m progress 
in describing the results of quench simulations in §]. 

To summarize, we have shown that the slow expansion of the matter in a relativis- 
tic nuclear collision forms disoriented chiral condensates through "annealing," rather than 
"quenching" Our methods are schematic and intended only to illustrate the underlying 
physics. Nevertheless, we argue that annealing can lead to domains significantly larger than 
those expected from a quench p. Our results therefore strengthen the suggestion that 
disoriented chiral condensates can be observed in heavy ion collisions at RHIC. 
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part of this work and to Gordon Baym and Miklos Gyulassy for discussions. B.M. thanks 
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work was supported in part by the Department of Energy grants DI5-AC02-76CH00016 and 
DE-FG05-90ER40592. 
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Figure Caption 



Figure 1: The Hartree effective potential as a function of $ = (ex, 7?) for the temperatures 
Tc = V^v and Tj3 compared to the T = potential Vq = A($2 _ ^2^4 

Figure 2: The behavior of the field for initial conditions (a) 7rO(r,) = 48 MeV, d7i^/dr{r,) = 
95 MeV/fm and (b) 71^{tc) = 5 MeV, dn^/dT{Tc) = 0. The time at which the system first 
drops below Tc is Tc = 10 fm. Dark curves are calculated for e = 1 and fight curves for 
6 = 0.5. 
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